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Abstract 

Stochastic invariant manifolds are crucial in modelling the dynam- 
ical behavior of dynamical systems under uncertainty. Under the 
assumption of exponential trichotomy, existence and smoothness of 
center manifolds for a class of stochastic evolution equations with lin- 
early multiplicative noise are proved. The exponential attraction and 
approximation to center manifolds are also discussed. 

Mathematics Subject Classifications (2010) Primary 60H15, 35R60; 
Secondary 37D10, 34D35. 

Keywords Stochastic partial differential equations (spdes), exponential 
trichotomy, center manifolds, stability, dynamical approximation. 

Contents 

1 Introduction 2 

2 Exponential trichotomy 5 

3 Stochastic center manifolds 8 



*Beijing International Center for Mathematical Research, Peking University, Beijing, 
China, mailto : chenxiao002214336@yahoo . cn 

^School of Mathematical Sciences, University of Adelaide, Adelaide, Australia. 
mailto : anthony . robertsOadelaide . edu . an 

•l-Institute for Pure and Applied Mathematics (IPAM), University of California, Los 
Angeles, USA & Department of Applied Mathematics, Illinois Institute of Technology, 
Chicago, IL 60616, USA. mailto:jduan@ipajn.ucla.edu 



1 



4 



Smoothness of center manifolds 



13 



5 



Exponential attraction principle 



20 



6 



Approximation of the center manifolds 



23 



7 



Examples show our results 

7.1 Stochastic reaction-diffusion modelling 

7.2 Semilinear damped wave dynamics . . 

7.3 An infinite dimensional center manifold 



26 

26 
28 
29 



1 Introduction 

Invariant manifolds are some of the most important invariant sets in non- 
linear dynamical systems. Stable, unstable, center and inertial manifolds 
have been widely studied in deterministic systems. But in many applications 
the nonlinear djTiamical systems are influenced by noise. Recently invariant 
manifolds for stochastic differential equations and stochastic partial differen- 
tial equations have been explored. Mohammed and Scheutzow [24] focused 
on the existence of local stable and unstable manifolds for stochastic differ- 
ential equations driven by semimartingales. Boxler [7] obtained a stochastic 
version of center manifold theorems for finite dimensional random dynam- 
ical systems by using the multiplicative ergodic theorem. Arnold [2] sum- 
marised various invariant manifolds on finite dimensional random dynamical 
systems. Roberts [30] assumed existence of stochastic slow manifolds for in- 
finite stochastic partial differential equations in explorring the interactions 
of microscale noises and their macroscale modelling. An unresolved prob- 
lem is the existence and nature of stochastic invariant manifolds on infinite 
dimensional spaces. 

Recently there are some theory developments about the problem. Duan 
and others [17, 18] presented stable and unstable invariant manifolds for a 
class of stochastic partial differential equations under the assumption of ex- 
ponential dichotomy or pseudo exponential dichotomy. Inertial manifolds 
are generalization of center-unstable manifolds on finite dimensional spaces 
to infinite dimensional spaces. Stochastic inertial manifolds on infinite di- 
mensional spaces are constructed by different methods [5, 11, 27]. Center 
manifolds are important invariant manifolds. It has been investigated exten- 
sively for infinite dimensional deterministic systems [4, 9, 34]. But there are 
only few papers dealing specifically with infinite dimensional center mani- 
folds. Gallay [20] proved the existence of infinite dimensional center-stable 
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manifolds and local center manifolds for a class of deterministic evolution 
equations in Banach spaces. 

This article establishes existence and properties of stochastic center man- 
ifolds on infinite dimensional Hilbert spaces with multiplicative white noise 
in the following class of stochastic evolution equation 

du 

— = Au + F{u) + uoW{t), (1) 

where m G if is a Hilbert space typically defined on some spatial domain, 
W{t) is the standard M- valued Wiener process on a probability space J^, P), 
which is only dependent on time. A simple example is the linear stochastic 
parabolic equation 

Ut = u^,: + Au + uoW{t), u{0,t) = u{Tr,t) = 0, t>0. (2) 

In this example the Hilbert space H is the function space L'^{0, it). This ex- 
ample stochastic partial differential equation (spde) has all three subspaces 
of interest: stable, center and unstable subspaces. Thus, the theory of Duan 
et al. [17, 18] and Caraballo et al. [12] does not apply. Sections 2-4 estab- 
lish the theory of center manifolds of such spdes. That is we prove the 
existence and smoothness of the center manifolds for stochastic evolution 
equation (1). The results generalize the work of Duan et al. [17, 18], who 
proved the existence and smoothness of center-stable and center-unstable 
manifolds. The stochastic center manifolds we obtain are generally infinite 
dimensional. Then Section 5 and 6 proves the principle of exponential at- 
traction from the stochastic center manifolds and an approximation to the 
stochastic center manifolds. 

Our main assumption is that there exists a spectral gap for the linear 
operator A. Such a gap allows us to construct the infinite dimensional center 
manifolds. Exponential dichotomy is one of the basic assumptions in non- 
linear dynamical systems. The assumption plays a central role in the study 
of stable and unstable manifolds. But the assumption of an exponential di- 
chotomy [17, 18] is not a sufficient condition for proving the existence of center 
manifolds. For example, a generator of a strongly continuous semigroup has 
some of its spectrum on the imaginary axis, which implies the existence of 
center manifolds, but Priiss [28] showed there is no exponential dichotomy. 
The pseudo-exponential dichotomy [17, 26] is also not a sufficient assump- 
tion to prove the existence of center manifolds since in the finite dimension 
case we need to ensure the existence of zero real part of the spectrum of a 
bounded operator. The concept of exponential trichotomy is important for 
center manifold theory in infinite dimensional dynamical systems and non- 
autonomous systems [3, 20, 26]. The existence of exponential trichotomy 
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means that the space is spht into three subspaces: center subspace, unstable 
subspace and stable subspace. In finite dimensional space the exponential 
trichotomy is satisfied automatically. Center manifold theory is based on the 
assumption of exponential trichotomy [26]. Section 2 introduces the expo- 
nential trichotomy and conjugated random evolution equations. Using the 
Lyapunov-Perron method, Section 3 proves the existence of stochastic center 
manifolds under the assumption of exponential trichotomy. We consider the 
existence of local center manifolds via a cut-off technique. When the stable 
subspace or the unstable subspace disappear, then our theorems reduce to 
the results by Duan and others [17, 18]. 

Deterministic center manifold theory is widely studied and further ex- 
tended by others in various contexts [1, 4, 9, 14, 21, 22, 34, 35, e.g.]. The 
method of proving the existence of center manifolds is to show that the center 
manifold is the fixed point of a certain of operator. It is also the intersection 
of the center-unstable manifolds and center-stable manifolds [1, 15]. In this 
paper, we define the stochastic center manifold as a graph of Lipschitz map. 
We directly prove the existence of such map by contraction mapping theo- 
rem, which is the method of Lyapunov-Perron [35]. It is different from the 
Hadamard's graph transform method [7, 18], which is more directly geomet- 
rical. Section 4 proves the center manifolds are smooth by using the method 
of Lyapunov-Perron. 

There are many applications of the theory of center manifolds such as the 
dimensional reduction [19, 37], bifurcation [9], discretisation [29, 31]. The 
stability of an equilibrium point on center manifolds has been considered 
in deterministic case [1, 25]. In the random case we show the equilibrium 
points are asymptotically stable when the equilibrium points restricted on 
the stochastic center manifolds are asymptotically stable (Section 5). We 
then prove that there exists an approximation of the stochastic centre man- 
ifold (Section 6). There has been some recent results on the approximation 
of invariant manifolds. Wang and Duan [3(:)] proved an asymptotic complete- 
ness property for a class of stochastic partial differential equations. Sun et 
al. [33] showed that an invariant manifold is approximated by a determin- 
istic manifold when the noise is small. Chen et al. [13] gave a geometric 
shape of invariant manifolds for a class of stochastic partial differential equa- 
tions. Blomker and Haire [6] considered a different approach to approximate 
stochastic partial differential equations by amplitude equations. Section 6 
gives a class of stochastic approximation to random evolution equations by 
stochastic center manifolds, which make some progress in the computation 
of stochastic center manifolds [30]. For this pathwise approximation version, 
the difficulty is to find a random variable to approximate the center mani- 
folds. This is different from the random norm [7] when using Carr's result [9]. 
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Section 7 uses some examples to illustrate our results. 

2 Exponential trichotomy 

We consider the stochastic evolution equation (1). We assume that the linear 
operator A : D{A) — )■ H generates a strongly continuous semigroup S{t) := 
exp{At) on a Hilbert space [H, | ■ |), which satisfies the exponential trichotomy 
with exponents a > 7 > > —7 > —/3 and bound K. 

First we introduce the exponential trichotomy, which generalizes the ex- 
ponential dichotomy [16, 17, 18, 32]. 

Definition 1. A strongly continuous semigroup S(t) := exp{At) on H is said 
to satisfy the exponential trichotomy with exponents a > 7 > > —7 > —f3 , 
if there exist continuous projection operators P'^, and P'^ on H such that 
the following four conditions hold. 

(i) id = P'^ + P'^ + P", where P'P^ = for i 7^ j, i,j G {c, s, u}, and id is 
the identity operator. 

(ii) P^exp(At) = exp(At)P^ P"exp(At) = exp(At)P", P"exp(At) = 
exp(y4t)P" fort>0. 

(Hi) Denote the reducing subspaces H"^ := P'^H , := P^H and iP* := 
P^H . We call H'^, and H'^ the center suhspace, the unstable sub- 
space and the stable subspace, respectively. The restriction exp{At)\Hi , 
t >0 are isomorphism from W onto H\ and we define exp[At)\Hi for 
t < as the inverse map of exp{—At)\Hi , i G {c,s,u}. 

(iv) 

\ex^{At)P%\<Kex^{-i\t\)\P%l teR,veH, (3) 
\exp{At)P''v\ < Kexp{at)\P''v\, t<0,veH, (4) 
|exp(At)P"t;| < i\:exp(-/3t)|P"^;|, t>0,veH. (5) 

Remark 1. If if is a finite dimensional space and there exist eigenvalues 
with real part zero, less than zero, and greater than zero, then we have 
the exponential trichotomy. If H is an infinite dimensional space and the 
spectrum of A satisfies (j{A) = o"^ U cr"^ U cr", where o"' = {A G cr{A) : Re A < 
a'^ = {A G a{A) : |ReA| < 7}, ct" = {A G a{A) : ReA > a}, and 
A generates a strong continuous semigroup, then we have the exponential 
trichotomy [20]. 
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We assume the nonlinear term F satisfies -F(O) = , and assume it to be 
Lipschitz continuous on iJ, that is, 

\F{ui)- F{u2)\ < LipFlm -n2| 

with the sufficiently small Lipschitz constant Lip F > . If the nonlinear 
term F is locally Lipschitz, let F^^\u) = xfi{u)F{u), where Xr(^) is a cut-off 
function, and then F^^^ is globally Lipschitz with Lipschitz constant i?LipF. 

Second we use a coordinate transform to convert the stochastic partial 
differential equation (1) into a random evolution equation [10, 17, 18]. Con- 
sider the linear stochastic differential equation, 

dz + fizdt = dW, (6) 

where /i is a positive parameter. A solution of this equation is called an 
Ornstein-Uhlenbeck process. Let Co(l^; 1^) be continuous functions on M, the 
associated distribution P is a Wiener measure defined on the Borel-a-algebra 
B{Co{R, M)). Define {Ot}tm to be the metric dynamical system generated by 
the Wiener process W{t). A random dynamical system on a metric space H 
on (Q, B, P, 6t) is a measurable map 

(f):Rx H xn H 

{t,u,u) (-7- (f)(t,U,u) 

such that 

(i) 0(0, = id; 

(ii) (f){t+s, u, u) = (f){t, u, 6sU})(j){s, u, u), for all t, s G M for almost all u & Q. 

Caraballo et al. [10] and Duan et al. [17] established the following lemma, 
which is used in the proof of the existence of stochastic invariant manifolds. 
Lemma 2. (i) There exists a {9t\t<m-'i'™0'™i^i set Vl G i3(Co(M, M)) of 
full measure with sublinear growth: 

lim M^ = o, CO en. 

i— >-±oo t 

(ii) For u G the random variable z(u) = —fi exp(/ir)a;(r) (ir exists 
and generates a unique stationary solution of the Ornstein-Uhlenbeck 
SDE (6) given by the convolutions 

.0 

Q X R 3 {u,t) ^ z{9tu) = —fi exp{fiT)9tu{T) dr 

J —oo 
.0 

= —fi / exp{fiT)uj{T + t) dr + fiuj{t). 

J — oo 

The mapping t (-)■ z{6tuj) is continuous. 
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lim l<^ = o forujen. 



(Hi) In particular, 

t— >±oo 

(iv) In addition, 

1 /■* 

lim - / z^Ot-u) dr = for u ^ Q . 

t— )>±oo t Jq 

We now replace the Borel-a-algebra S(Co(M, M)) by 

J^={nnD:De B{Co{R, R))} 

for Q given in Lemma 2. The probability measure is the restriction of the 
Wiener measure to this new cx-algebra, which is also denoted by P. In the 
following we consider a random dynamical system over the metric dynamical 
system J^,P,e). 

We show that the solution of equation (1) define a random dynamical 
system. For any u* E H and a; G , we introduce the coordinate transform 

u = T{u*, uj) = u* eyc^lziOtOj)] (7) 

and its inverse transform 

T^^{u,uj) = uexp[—z{9tuj)]. (8) 

Without loss of generality, we drop the *. We transform the stochastic 
evolution equation (1) into the following partial differential equation with 
random coefficients 

du 

— = Au + z{9tu;)u + G{9tu;,u), m(0) = % G if, (9) 

(JjL 

where G{u,u) = exp[—z{u)]F{exp[z{uj)]u). We have 

\G{u,ui) - G{u,U2)\ < Lip„G \ui - M2I, 

where Lip„ G denotes the Lipschitz constant of G{-, u) with respect to u. For 
any u E VL the function G has the same global Lipschitz constant as F by 
the construction of G. Define 



'^A{t,s) = exp 



A{t- s)+ / ziOru) dr 
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as a 'state transition operator' for the linear random partial differential equa- 
tion (9) in the case when the nonlinearity is neglected, G{u, u) = . The 
solution is interpreted in a mild sense 

u{t) = •^A{t,0)uo+ I <l/A{t,s)G{6sUJ,u{s))ds . 
Jo 

This equation has a unique measurable solution from the Lipschitz conti- 
nuity of G and measurable property of z{u). Hence the solution mapping 
(t,v,u) H-> u(t,v,uj) generates a random dynamical system. 

Duan et al. [18] proved that the solution of (9) is the coordinate trans- 
formation the solution of (1). 

Lemma 3. Suppose that u is the random dynamical system generated by (9). 
Then 

{t, V, u) I— 7- T~^{9tUj, u(t, T{uj, v),uj)) =: u(t, v, u) (10) 

is a random dynamical system. For any v & H this process {t, u) h-> u{t, v, u) 
is a solution to (1). 



3 Stochastic center manifolds 

We introduce the definition of stochastic center manifolds. A basic tool of 
proving the existence of stochastic center manifolds is to define an appropri- 
ate function space, which is a Banach space. 

Definition 4. A random set M{uj) is called an (forward) invariant set for a 
random dynamical system (j){t,u,x) if (f){t , M (u) , u) C M{6tu) for t > 0. 
If we can represent Miuj) by a graph of a (Lipschitz) mapping from the 
center subspace to its complement, h'^{-,uj) : H'^ — )■ © , such that 
M{uj) = {v + h'^{v,uj) : V G H'^}, h'^{0,uj) = , and the tangency condition 
that the derivative Dh'^{0,uj) = 0, h'^{v,-) is measurable for every v G H'^, 
then M{uj) is called a (Lipschitz) center manifold, often denoted as M'^(co'). 
M{uj) is called a local center manifold if it is a graph of a (Lipschitz) map- 
ping XR{.v)h''{-,uj). 

We first show the existence of a Lipschitz center manifold for the random 
partial differential equation (9). Then, we apply the inverse transforma- 
tion to get a center manifold for the stochastic evolution equation (1). 

For each 77 > , we denote the Banach space 



Gr, 



H) : supexp 



?7|t| — / z(9rCo) dr 



m)\<oo 
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with the norm 



sup exp 



-ri\t\ 



\m- 



The set Cr, is the set of 'slowly varying' functions. We know that the functions 



are controlled by exp 



?7|t| + Jg z{9rOj) dr 



Let 



M''{uj) = {uq E H : u{-,Uo, tu) E C,,} , 

where u(t, uq, u) is the solution of (9) with the initial data u{0) = uq . 

We prove that M'^{uj) is invariant and is the graph of a Lipschitz function. 
For simple the expression of the proof, define the operator B : M. L{H) 

by 

'-exp(At)P", t<0; 
exp(yl^)P^ t > , 



B{t) 



Different from the proof process of Duan et al. [17, 18], we need analysis the 
behavior of the solution on center subspace. 

Theorem 5. If < rj < min{/3,a} such that the nonlinearity term is suffi- 
ciently small, 



77 — 7 (3 — rj a — 7] 



<1, 



:ii: 



then there exists a center manifold for the random partial differential equa- 
tion (9), which is written as the graph 

M\uj) = {v + h'{v,uj) : V E iP}, 

where h'^{-,u) : H'^ — )■ if" © is a Lipschitz continuous mapping from the 
center subspace and satisfies h'^{0,uj) = . 

Proof. First we claim that Uq E if and only if there exists a slowly 

varying function u{-,Uo,uj) E Cn with 



u{t,uo,uj) = '^A{t,0)v [ '^A{t,s)P''G{6sUJ,u{s))ds 

Jo 

/+00 
^o(^, s)B{t - s)G{esUj, u{s)) ds , 
-00 

where v = P'^uq, \I/o(^, s) = exp z{9rOj) dr . 



(12) 
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To prove this claim, first we let Uq G M'^{u). By using the variation of 
constants formula, the solution on each subspace denoted as 

P^u(t,uo,u) = '^A{t,0)v + [ '^A{t,s)P^G{esU,u)ds . (13) 

Jo 

P^'uit, uo, u) = ^A(t, r)P"n(r, uq, ^) + ^ ^^(t, s)P"G(e,a;, u) ds . (14) 

P'u{t, Uo, u) = ^A{t, t)P'u{t, no, ^) + ^ ^A{t, s)P'G{e,uj, u) ds . (15) 

Since the slowly varying function u G C,; , we have for t < r that the magni- 
tude 

\^A{t, t)P''u{t, Uo, uj)\<K exp[a{t - r)]^oit, 0) exp{r]T)\u\c^ 

= K^ait, 0) exp[-(a - v)^]\u\c^ 
— >■ as r — +oo. 

For t > T , 

\<i/A{t,r)P'u{T,Uo,uj)\ < K exp[-(3{t - T)]^o{t,0) exp{7]T)\u\c, 

= K^_^{t,0)exp[{l3 + 7])T]\u\c^ 
— >■ as r — )■ — oo. 

Then, taking the two separate limits r — )■ ±oo in (14) and (15) respectively, 
P"u(t,uo,a;) = / ^A{t,s)P''G{esUj,u{s))ds, (16) 

J oo 

P'u{t,uo,uj)= I ^A{t,s)P'G{esUj,u{s))ds. (17) 

J — oo 

Combining (13), (16) and (17), we have (12). The converse follows from a 
direct computation. 

Next we prove that for any given v G H'^, the centre subspace, the integral 
equation (12) has a unique solution in the slowly varying functions space C^. 
Let 

J\u,v):=^A{t,Q)v+ I ^A{t,s)P^G{esUJ,u{s))ds 

(18) 



+ 00 



^o(t, s)B{t - s)G{esU, u{s)) ds . 
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is well-defined from x H'^ to the slowly varying functions space C^. For 
each pair of slowly varying functions u,u E C,^ , we have that for 7 < 77 < 
min{/3, a}, 

\r{u,v)-.r{u,v)\c. 



< sup < exp 



-?7|t| 



ziOsUj) ds 



^^(t, s)P\G{esUJ, u) - GiOsUJ, u)) ds 



+ 



/+00 
*o(t, s)B{t - s){G{esUj, u) - GiOsio, u)) ds 
-00 

/ exp[{-f - ri)\t - s\] ds 
Jo 



< sup < KLip^ G\u — u\ 



+ / exp[(r/-/3)(t-s)]rfs + 
1 1 1 



/+00 1 
exp[{a — ri){t — s)] ds > 



<i^Lip„G 



77 — 7 (3 — 1] a — 7] 



1'" - ^Ic,- 



(19) 



From equation (19), J'^ is Lipschitz continuous in v. By the theorem's pre- 
condition (11), J'^ is a uniform contraction with respect to the parameter v. 
By the uniform contraction mapping principle, for each v G H"^, the mapping 
J'^{-,v) has a unique fixed point u{-,v,uj) G . Combining equation (18) 
and equation (19), 

K 



\u{-,v,uj) - u{-,v,u)\c^ < 



1-KUp^g(^ + ^ + ^ 

\ P^V a—Tj 



\v — v\ 



(20) 



for each fixed point m(-, f , u). Then for each time t > , u{t, ■, u) is Lipschitz 
from the center subspace to slowly varying functions C^. u{-,v,u) G C,, 
is a unique solution of the integral equation (12). Since u{-,v,u) can be an 
tu-wise limit of the iteration of contraction mapping J'^ starting at and J'^ 
maps a ^-measurable function to a J-'-measurable function, u(-,w,a;) is J-"- 
measurable. Combining u{-,v,u) is continuous with respect to if, we have 
m(-,w,cj) is measurable with respect to (■,f,a;). 
Let h^{v,uj) := P"m(0,w,u;) ©P"M(0,w,a;). Then 

r" + 00 



^0(0, s)B{-s)G{6sU}, u{s, V, u)) ds . 



We see that h'^ is J-'-measurable and h'^{0,u) = 0. From the definition 
ofh'^{v,uj) and the claim that mq G (cu ) if and only if there exists u{-,uq,u) G Gr, 
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with n(0) = uq and satisfies (12) it follows that uq G M'^{u) if and only if 
there exists v G H'^ such that uq = v + h'^{v,u), therefore, 

M^{uj) = {v + h%v,u) : V G H'}. 

Next we prove that for any x & H , the function 

inf \x - {y + h''{y,u))\ (21) 

is measurable. Let H' be a countable dense subset of the separable space H. 
From the continuity of h'^{-,u), 

inf \x-{y + h''{y,uj))\= inf \x - P^y - h'iP'y^uj)]. (22) 

The measurability of (21) follows since u — h'^{P'^y, u) is measurable for any 
yeH'. 

Finally, we show that M'^{u) is invariant, that is for each uq G M'^{u), 
u{s,uq,u) G ]iP^{9sUj) for all s > . Since for s > , u{t + s,Uo,u) is a 
solution of 

du 

— = Au + z{6t{6sUj))u + G{6t{6su), u), n(0) = u{s, uq, u). 

Thus M(t, m(s. Mo, w), 6^30;) = u(t + s, mq, w) and u{t, u{s, uq, u), OgUj) ^ Cr, ■ So 
we have u{s,uo,u) G M'^{9sUj) 

□ 

Theorem 6. Suppose that u is the solution of random partial differential 
equation (9), then M^{u) = T~^{u, M'^{uj)) is a center manifold of the 
stochastic partial differential equation (1). 

Proof. Denote u{t,uj,x) the solution of (1). From Lemma 3, for t > , 

u{t, u,M\uj)) = T^\etUJ, u{t, 00, T{uj, M'(w)))) 

= T-\etUj,u{t,uj,M'{uj))) (lT-\etUJ,M\etUj)) =M\etUj). 

So M'^{uj) is an invariant set. Note that 

M\u) = T~\u,M\uj)) 

= {no = T-\uj, V + h%v, uj):ve H^] 

= {no = e"(^)(M + /i"(M,w)) : v G H'} 

= {no = M + e"('^)/i"(e-"('^)M,w) : v G i/"}, 

which implies that M'^{u) is a center manifold. □ 
Note that a local center manifolds is not unique [9]. 
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4 Smoothness of center manifolds 

In this section, we prove that for each u E Q , M'^{u) is a center manifold. 
Theorem 7. Assume that G is in u. If < krj < min{/3,a} and 

KLip^Gi — — + -7^ — + — ^ — )<1 foralll<i<k, 
\iri — 7 p — iTj a — If] ) 

then M'^{uS) is a G^ center manifold for the random evolutionary equation (9), 
and D/i^(0,w) = . 

Proof. We prove this theorem by induction. First, we consider k = 1 . Since 

i^Lip„Gf^ + -^ + ^— ) < 1 

\?7 — 7 p — rj a — rj J 

there exists a small number S > such that < t] — rj' < min{/3, a} and for 
all < r/' < 25 , 



KLip^G 



+ -P. -, ^ + 



{rj — rj') — (3 — [rj — rj') a — {rj — rj') 



< 1. 



Thus, v) defined in the proof of Theorem 5 is a uniform contraction in 
C^_^' C Grj for any < rj' < 26 . Therefore, u{-, v,u) G C^_,,' . For vq e H'^, 
we define two operators: let 

Svo = ^A(t,0)Wo, 

and for v G Grj^s let 
ft 



Tv = f •^A{t,s)P''DuG{6sUJ,u{s,vo,Lo))vds 
Jo 



+ 00 



+ / "^olt, s)B{t — s)DuG{6sUJ,u{s,vo,u))v ds . 

From the assumption, is a bounded linear operator from center subspace H'^ 
to slowly varying functions space Gr/s- Using the same arguments in the 
proof Theorem 5 that J'^ is a contraction, we have that T is a bounded linear 
operator from C^_5 to itself and 

\\t\\<kup^g( — I — + — 7T + — r — k1 <1' 

V^-'^-7 /5-(r/-<5) a-{r]-6)J 
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which imphes that the operator id — T is invertible in C^s- For ^! £ 
we set 



— DuG{9sU!, u{s, Vq, u)){u{s, V, u) — u{s, Vq, u)) ds 



+00 



\l/o(t, s)B{t — s) G{9sUJ, u{s, V, u)) — G{9sUJ, u(s, vq, u)) 



- D^G{9sUJ,u{s,Vo,uj)){u{s,v,uj) - u{s,vo,uj)) 



ds 



We obtain 

U{-, V, U) - U{-, Vq, u) - T{u{-, V, u) - u{- , Vq, u)) = S{v - Vq) + I, (23) 

which yields 

V, u) — Vq, u) = (id —T)~^S{v — Vq) + (id —T)~^I. 

If |/|c^_^ = o{\v — t>o|) as f — )■ Vq, then u{-,v,u) is differentiable in v and 
its derivative satisfies Dyu(t,v,u) G L{H'^ , G^^^s) , where L{H^,G^^s) is the 
usual space of bounded linear operators and 

Dyu{t, V, cu) = \&a(^, Q)v + / '^Ai't, s)P''DuG{9sUj, u{s, v, uj))D^u{s, v, to) ds 



+00 



\l/o(t, s)B(t — s)DuG{9sUJ, u{s, v, u))D^u{s, v, u) ds . 

(24) 



The tangency condition Dh'^{0,uj) = is from equation (24). 
Now we prove that 

\I\c,,s = oi\v-vo\) 



(25) 



as f — )■ Vq . We divide I into several sufficient small parts. Let be a large 
positive number to be chosen later and define the following ten integrals 



Ii = exp 



iv-S)\t\ - / z{9sCo)ds 



N 



^A{t,s)P'' G{9sC0,u{s,v,u)) 



G{9sUJ,u{s,Vo,u)) - DuG{9sUJ,u{s,Vo,u)){u{s,v,u) - u{s,Vo,U!)) 



ds 



14 



for t > iV . 



I[ = exp 



-{V-S)\t\ - / z{9,co)ds 



ds 



for t< -N . 



I2 = exp 



iv-S)\t\ - / z{9sCo)ds 



N 



G{9sUj, u{s, vo, u)) - DuG{9sU}, u{s, vq, u)){u{s, v, u) - u{s, vq, to)) 



ds 



for < t < AT . 



I2 = exp 



-iv-W - / z{9sio)ds 
'0 



N 



- G{9sUJ, U{S, Vq, u)) - DuG{9sUJ, U{S, Vq, Uj)){u{s, V, u) - u{s, Vq, Uj)) 



ds 



for -iV < t < . 

Let be a large positive number to be chosen later. For |t| < , we set 



I3 = exp 



-{V-S)\t\ - / z{9,co)ds 



N 



^A(t,s)]P" Gi9sU,u{s,v,co)) 



G{9sU, U{S, Vq, u)) - DuG{9sU, U{S, Vq, V, u) - u{s, Vq, u)) 



ds 



I3 = exp 



-{r]-6)\t\- / z{9,co)ds 



N 



^A{t,s)P' G{9,u,u{s,v,u)) 



G{9sU, U{S, Vq, u)) - DuG{9sU, U{S, Vq, u)){u{s, V, u) - u{s, Vq, u)) 



ds 



I4 = exp 



ir]-6)\t\- / z{9su)ds 



N 



[^^(t,s)]P" G{9,u,u{s,v,u)) 



G{9sU},u{s,vo,u})) - DuG{9sUJ,u{s,vo,u)){u{s,v,u) - u{s,vq,uj)) 



ds 



I'^ = exp 



-{ri-S)\t\- / z{9,co)ds 



^Ait,s)P' Gi9sLo,u{s,v,u)) 



G{9sUJ,u{s,vo,u})) - DuG{9sUJ,u{s,vo,u)){u{s,v,u) - u{s,vo,uj)) 



ds 
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For |t| > iV, we set 



h = exp 



iv-S)\t\ - / z{9sCo)ds 



+00 



[^A{t,s)]P^ G{9sCo,u{s,v,uj)) 



G{dsUJ,u{s,vo,u)) - DuG{6sU},u{s,vo,uj)){u{s,v,u) - u{s,vo,uj)) 



ds 



5 = exp 



-{r]-6)\t\ - / z{esuj)ds 
Jo 



^Ait,s)P' G{9,tu,uis,v,tu)) 



- G{9sUJ,u{s,Vo,uj)) - DuG{9sUj,u{s,vo,uj)){u{s,v,uj) - u{s,Vo,uj)) 



ds 



It is sufficient to show that for any e > there is a cr > such that if 
|w — ^ol < cr 5 then |/|c^_^ < e\v — Vo\. Note that 

l-^lc _i ^ sup Ii + sup I2 + sup I[ + sup + sup Is + sup I4, + sup /s 

t>N N>t>0 t<~N -~N<t<0 it\<N \t\<N \t\>N 

+ sup I's + sup 14 + sup Jg. 

|t|<lV |i|<lV |t|>lV 

A computation similar to (20) imphes that 
Ji < 2Js:Lip„G' 

exp[(7 -{rj- 6))\t - s\] exp{~6\s\)\u{-,v,uj) - «(-, Wq, w)|c,_2i ds 
2K^ Lip„Gexp(-5Ar) 



N 
< 



(^-^-5)[l-KUp^G 
Choose so large that 



(r7-25)-7 ^ l3~{ri~2S) ^ a-{r]-25) 



V - Vo 



2K^ Lip„G'exp(-5iV) 



{7i-l-5)\l-KUp^G 



(r;-25)-7 ' ^-(r)-25) ' a-{ri-25) 

Hence for such N we have that 

sup/i < ^€\v -Vo\. 
t>N o 



1 

< -e. 
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Fixing such A^, for I2 we have that 
/2<ir^ exp[(7 + (r/-5))|t-s|]{^ \D^G{9sUJ,tu{s,v,co) + {1 - r) 



< 



u{s,Vo,uj)) - DuG{esUJ,u{s,Vo,uj))\dTj\u{-,v,uj) - u{-,Vo,uj)\c^_g ds 
K'^\v-vo\ 



1 - K Lip„ G 



rj—S—'y fi—{ri—&) ci—(rj—5) 

exp[(7 + (r/ - 5))\t -s\]{J^ \DuG{9sC0, tu{s, v, cu) + {1 - r) 

u{s, Vq, u)) — DuG{9sU, u{s, Vq, cu)) I dr| ds . 

From the continuity of the integrand in {s,v), the last integral is continuous 
at the point vq- Thus, we have that there is a ai > such that if |w — i^ol < cti , 
then ^ 

sup I2 < -e\v - vo\. 

N>t>0 o 

Therefore, ii \v — vo\ < ai , then 

1 , 

sup/i + sup I2 < -e\v - Vo\. 

t>N N>t>0 4 

In the same way, there is a cr^ > such that if if |f — WqI < then 

sup h + sup I2 < '^e\v - vo\. 

t<-N ~N<t<o 4 

Similarly, by choosing to be sufficiently large, 

sup_/4 + sup_/5 < ^e\v - Vol, 

\t\<N \t\>N ° 

sup_r^ + sui^/5 < ^e\v - Vol, 

\t\<N \t\>N ° 

and for fixed such A^, there exists a"2 > such that if |f — fo| < (^2 , then 
sup I3 < ^<^\vi - V2\ and sup Jg < ^e|fi - V2\- 

\t\<N ° \t\>N ° 

Taking a = minjcri, a[, (T2}, we have that if It" — fo| < o" , then 
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Therefore \I\c^_s = o{\v - vo\) as v ^ vq . 

We now prove that D^u{t, v, u) is continuous from the center subspace H'^ 
to slowly varying functions space C^. For v,Vo G H'^, using (24), 

Di,u{t, V, uj) — Di,u{t, t>o, u) 

^^(t, s)P'' {DuG{dsUJ, u{s, V, uj))D^u{s, v, u) 

— DuG{9sC0, u{s, Vq, u))D^u{s, Vq, Cj) j ds 

+ / ^Q{t, s)B{t — s)iDuG{9sU!,u{s,v,u))Dyu{s,v,u) 

J — oo 

- DuG{9sU, u{s, vo, uj))D^u{s, Vq, u)^ ds ^26) 
^A{t,s)P''(^DuG{esUj,u{s,v,u)) 

{Dyu{s, V, u) — Di,u{s, Vq, co)) ) ds 



+ I ^o{t,s)B{t - s)(DuGiesU,u{s,v,u)) 

J — oo 

{Dyu{s, V, cu) — Di,u{s, Vq, co))j ds + I, 

where 

1= [ ^A{t,s)P'{D^G{esU,u{s,v,u)) 
Jo 

— DuG{9sC0, u{s, Vq, Uj))Di,u{s, Vq, Co) ds 

/ + 00 
^o(t, s)B{t - s) {DuGie^uj, u{s, V, u)) 
-oo 

— DuGiOgUJ, u{s, Vq, Uj))DvU{s, Vq, u) ds . 

Then from estimating \Dyu{-,v,u) — Dyu{- , Vq, u)\l(^h<'^c^), we have 

mL{H'=,C^) 



\D^u{-,V,Uj) - DyU{-,Vo,Uj)\L(H-,C^) < 



i-kup^g(^^ 



, — + — 

v-j P-v <^~v , 



Using the same argument we used for the last claim of equation (25), we 
obtain that c^) = o{\v — vo\) as v vq . Hence Dyu(t,v,u) is contin- 

uous from the center space H'^ to bounded linear operators space L{H'^, Grj). 
Therefore, u{t,v,u) is G^ from H'^ to C^. 
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Now we show that u is C'^ from the center space H"^ to slowly varying 
functions space C^^ by induction for k > 2 . By the induction assump- 
tion, we know that u is C'^"^ from H'^ to C(^k-i)r] and the {k — l)st deriva- 
tive Dl-^u{t, v,u) satisfies the following equation. 



ct 

Dl^\= I '^A{t,s)P%nuG{9,u,u)D'^~^uds 
Jo 

/ + 00 
^o{t, s)B{t - s)D^G{esUJ, u)D^-\ ds 
-oo 

Jo 

/ + 00 
^oit, s)B{t - s)Rk-i{s, V, u) ds , 
-oo 



where 



1=0 ^ ^ 

Note that -D*n G Cj^ for z = 1, . . . , /c — 1 , from the induction hypothesis. 
Thus, using G is G^, we verify that Rk-i{-,v,uj) G L^~^[H'^,G(^k-i)r]) and 
is G^ with respect to v, where L^~^ {H^, C(fc_i)^) is the usual space of bounded 
k — 1 linear forms. Since 

JsTLip^G ( + — ^ — + — ^ — )<1 for all 1 < i < fc. 

\irj — 7 p — If] a — iTj J 

The same argument used in the case k = 1 shows that D'^~^u{-,v,uj) is G^ 
from H'^ to L^{H'^,Gkri)- This completes the proof. □ 

Theorem 8. Assume that the nonlinearity F{u) in equation (1) is G^ smooth. 
If ^ < krj < mm{P,a} and 

KLipFl — - — + — ^ — + — - — )<1 foralll<i<k, 
\iri — 'J p — irj a — ir] J 

then M'^{u) = T^^{uj, M'^{u)) is a G^ center manifold for the stochastic 
partial differential equation (1). 

Proof. From the construction of center manifolds, we obtain 

M\uj) = {v + h%v,u) -.veH'}, 

where h'^{v,u) = exp[z{u)]h'^{exp[—z{u)]v,uj). Note that h'^{-,u) is G'', then 
h'^{-,oj) is G^. We conclude M'^(a;) is a C^' center manifold for the stochastic 
partial differential equation (1). 

□ 
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5 Exponential attraction principle 

This section proves the stabihty of the random evolution equation (9). Our 
main result is Theorem 10, which shows that the dynamic behavior of equa- 
tion (9) is exponentially approximated by the solution of on stochastic center 
manifold. We do not assume H'^ = or iJ" = . Instead, we project the solu- 
tion of equation (9) to the stable and unstable subspaces. Lemma 9 is a week 
approach solutions version as it only concerns solutions approach M'^{uj). It 
does not assert solutions. 

Lemma 9. Let M'^{uj) be a locally stochastic center manifold for the stochas- 
tic evolution equation (9), then there exist positive constants U and fi such 
that for u ^ Q , 

\P'u{t, uo, uj) - P'h\P%{t, uo, uj),uj)\<U exp{-i2t)\P'uo - P'h\P^uo, u) \ , 
t > large enough. 

Proof. Let X{t, uq, oj) = P^u{t, uq, uj) — P'^h'^{P'^u(t, uq, uj),u), then X satis- 
fies the equation 

dX(t, Uq, u) 



dt 

where 



AX + z{etUj)X + N{X, Otuj) , (27) 



N{X, Otu) = P'G{9tU, (id -P')u + X + P'h%P^u{t, uj),uj)) 

- P'G{9tUj, (id -P')u + P'h'iP'uit, Mo, w), w)) 
+ DP'h%P%{t, no, uj)),u) [P'G{etUJ, (id -P')u 

+ P'h\P^u{t,UQ,Uj),Uj)) 

- P^G{etU, (id -P')u + X + P'h%P%{t, Mo, 

Since the derivative of stochastic center manifold is less than Lipschitz con- 
stant Lip^ G and h'^ is Lipschitz continuous, we have \N{uj,X)\ < Ui\X\ for 
a positive constant Ui. We obtain 

X{t,uo,u) = *^(t,0)X(0,Mo,w) + / ^A{t,s)N{X,u)ds . 

Jo 

The conclusion of Lemma 2 yields z{9ru)dr < et , z{9rUj)dr < e{t — s) 
for t large enough and /3 > e > . Then 

\X{t,uo,uj) \< Kexp[{-f3 + e)t]\X{0,uo,uj)\ 

+ UiK / exp[{-(3 + e){t - s)\X{s,uo,uj)\ds . 
Jo 

From Gronwall's inequality we have the lemma. □ 
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Remark 2. Similarly we have the conclusion there exist positive constants V 
and u such that for u ^ Q, 



\P''u{t,uo,uj) - P''h'{P'u{t,Uo,uj),cu)\ < 1/exp(z/t)|P"no-P"/i^(w,w)|, 

t < small enough. In the finite dimensional case [7, Theorem 7.1], the 
conclusions are expressed as the random norm. 

Now we consider the asymptotic behavior of the stochastic evolution equa- 
tion (9). The stochastic dynamical system on the stochastic center manifold 
is governed by the following random evolution equation. 

V = P^Av + z{etuj)v + P^G{etUJ, V + h\v, uj)). (28) 

If the subspace H'^ is finite dimension, then equation (28) is an ordinary dif- 
ferential equation. The following theorem shows that the solution of stochas- 
tic evolution equation (9) is exponentially approximated by the solution of 
equation (28). 

Theorem 10. Let M'^{u) be a locally stochastic center manifold for stochastic 
evolution equation (9), then there are positive random variables Liiuo) and L2{uj) 
and constants li and I2 such that for initial values Uq & H and Vq G M^, the 
following two conditions hold: 

\P'^u(t,uo,u) — v{t,vo,uj)\ < Li{u) exp{—lit) , (29) 
\P'u{t,uo,uj) - P'h%v{t,Vo,u),u)\ < 12(00) expi-kt), (30) 

t > large enough and u & Q. 

Proof. We construct a differential equation about the error of the two solu- 
tions on centre subspace. Let X'{t,uo,oj) = P'^u(t,uo,uj) — v(t,vo,u), which 
is the errors of the two solutions on centre subspace, then X' satisfies the 
equation 

dX'{t, uq, u) 
Jt 

where the nonlinear term 



AX' + z{etu)X' + N'{X',u), (31) 



N'{X', oj) = P'G{9tUJ, u{t, uo, to)) - P'GiOtU, v + h'iv, to)). 

Rewrite the nonlinearity A^' as three steps, so that A^' is bounded. 

N'{X', to) = P^G{6tU0, u{t, no, to)) - P^'GiOtto, P'u + /i'=(P%, to)) 

+ P'GiOtto, P'u + /i"(P^M, to)) - P^GiOtto, P^u + h\v, to)) 
+ P'G{etto, P'u + h'^iv, to)) - P^G{etto, V + h'iv, u)). 
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Rewrite equation (31) as an integral equation. 

X'{t,Uo,u) = ^A{t,0)X'{0,uo,u)+ [ ^A{t,s)N'{X',u)ds. 

Jo 

Since G is Lipschitz continuous, there exist positive constants Vi, vi and V2 
such that 

\N'{X',u)\ <Viexp{-vit){\P'uo - P'h''{P''uo,uj)\ + |P"mo - ^"/^'(P'mo, w)|) 

+ V2\X'\. 

Lemma 2 yields z{9rUj)dr < et and z{9rUj)dr < e{t — s) for t large enough 
and 7 > e > . It follows that 

\X'{t,uo,co)\ <Kexp[{-f + e)t]\X'{0,uo,co)\ 

-K{^ + e + vi)-^ exp{-vit)Vi{\P'uo - P'h'{P'uo, u)\ 

+ V2K [ exp[{-f + e){t-s)\X'{s,Uo,uj)\ds. 
Jo 

From Gronwall's inequality we have inequality (29). 
Now we prove inequality (30), 

\P''u{t,Uo,u) - P''h''{v{t,vo,u),u)\ 

< \P'u{t,uo,u) - P'h\P''u,u)\ + \P'h\P^u,u) - P'h\v{t,vo,u),u)\. 

By Lemma 9 and inequality (29) we have inequality (30). □ 

Remark 3. Similarly we conclude there exist positive random variable ^3(0;) 
and positive constant /a such that 

|P'^M(t,Mo,u;) - ^7(^,^70,^^)1 < ^3(w)exp(/3t), (32) 
\P''u{t,UQ,uj) - P''h\v{t,VQ,uj),uj)\ < L^{uj)exp{ht), (33) 

t < small enough and a; G fi. 

Corollary 11. Let = 0. Suppose that the zero solution of equation (28) 
is asymptotically stable and the initial value Uq is sufficiently small. Then 
the zero solution of equation (9) is asymptotically stable. 

Proof. From Theorem 10, we conclude that the zero solution of equation (28) 
is asymptotically stable on centre and stable subspaces respectively. This 
establishes Corollary 11. □ 
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6 Approximation of the center manifolds 

In this section, we prove the existence and approximation of stochastic center 
manifolds. We avoid assuming the noise is different iable in the proof pro- 
cess [8]. Denote h{v{t,Vo,u),u) =: cu) for the initial value P'^M(t, no, w) =: 
v{t, Vq, u). For a J-'-measurable function g{v, u) : H'^ — H^(BH^ which are 
in a neighbourhood of origin in H'^, define 

M{g){v, u) = gt{v, to) + g^,{v, lo) [P^'Av + z{etio)v + P^G{9tUj, v + g)] 
-Ag- z{etuj)g - P'G{etCO, v + g) - P"G(^ia;, v + g), 

where v G H'^. 

Theorem 12. Suppose that g is a J-" -measurable function such that g{0,uj) = 
, gv{0, u) = and M{g){v, u) = |'^) for some q > 1 . Then as v , 

\h''iv,u)-giv,u)\ = 0{\v\''). 

Proof. Let an J^-measurable function y{s,v,u) G C^(w) and let y{s,-,u) : 
TP^ — © H'^ be continuous differential with compact support such that 
y{0,0,u) = g{0,u) = 0, yy{0,0,u!) = gv{0,uj) = and yy{0,v,uj) = gv{v,u) 
for |f| small enough. The function y{s,v,u) connects the functions h'^{v,u) 
and g{v,u). Set 

M{y){0,v,io) = yt{0,v,io) + y,,{0,v,io)[P'Av + z{etuj)v 

+ P'GiOtUJ, y{0, V, u))] - Ay{Q, v, oo) - z{etuj)y{0, v, u) 
- P'GiOtu, 2/(0, V, u)) - P'^GiOtu, y(0, v, u)), 

where v G TP^. For a J-'-measurable function w{s,v,uj), define the operator 

Sw = J^{w{s, V, u) + y{s, V, u!),v) — y{s, v, u), 

where the operator J'^ is defined by 

J''{u,v) := exp{As)v+ f exp[A{s - r)]P''G{6rU),u{r)) dr 

Jo 



+ 00 



« (34) 
B{s - r) exp( / z{9iuj) dl)G{erUj, u{r)) dr . 

J r 



For a random variable K{u), define the domain Y 

Y{u) = {w{s,v,uj) G C^{uj) : \w{0,v,uj) \ < K{uj)\v\'' for all v e H""} . 

From its construction, Y{u) is closed in C^(a;). Since J'^ is a contraction 
mapping on C^(a;), 5 is a contracting mapping on Y{u). We prove that 
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there exists a random variable K{u) such that 5* maps Y{u) into Y{u). 
Then by the uniqueness of fixed points and definition of h'^ we conclude the 
theorem. 

Now we prove Y{u) exists. For w G Cr,{co), let Uc{s,v,u) be the solution 

of 

du 

— ^ = Auc + z{9sUj)uc + G{6sUj, y{0, Uc, u) + w{0, Uc, u)), 
at 

Uc{0, V, u) = V. (35) 
From equation (34) and the construction of h{v,u), 

5(-r)exp / z{9iu)dl 

■oo \_J r 

GiOr-u, y{0, Uc, uj) + w(0, Mc, w)) dr . 
Since the slowly varying function y{r,Uc,U!) G C^(ci;) , 



exp 



-Ar+ I z{9iu)dl 

J r 



< ii'exp[(/3 - ?7)r]|y(r,Mc,w)|c„(a;) 
— )■ as r — )■ — oo. 



Also 



exp 



Ar+ z{0iuj)dl 



< i^exp[-(a - r])r]\y{r, Uc,w)|c^(<^) 
— )■ as r — )■ +00. 



So 



y{0,v,u 
'•° d_ 
dr 

d 
dr 











exp 



exp 
exp 
-Ar 



Ar 



Ar 





ziOiUj) dl 



z{6iuj) dl 
z{6i(jj) dl 



P'y{0,u„u) 



dr 



dr 



+ -^P'y{0,u,,u) 
dr 



dr 



exp 



_ z{e.rUj))PMO, n„ oo) + — P"y(0, n„ u) 

dr 



A- z{erio))P'y{0,u,,uj) 



-Ar+ / z{0,uj)dl 



dr , 
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Here 

{-A - z{erUj))P'y{0, u,(r, v, u),u) + -^P'yiO, n,(r, v, u),u) 

ar 

= {-A - z{eru))P'y{0, Ucir, v, u),u) + P'VviO, u^r, v, u),u)[Auc 

+ z{9rUj)uc + P^G{6rUj, y{0, Uc, uj) + w{0, Uc, uj))] 

+ P^yr{0,Uc{r,v,uj),uj) 
= M{P'y) - P'y^{0,Ucir,v,uj),uj)[Auc + z{9ruj)uc 

+ P^GiOrU, 2/(0, u„ u))] + P^GidrU, P'y{0, u„ u)) 

+ P'yy{0,Uc{r,v,uj),uj) 

[AUc + z{9rUj)Uc + P''G{drUJ, 1/(0, w) + w(0, Oo))] 

= M{P'y) - P'y,{0, u,ir, v, u),u) [P^G(^,u;, y(0, u,, u)) 

- P^'GiOrUJ, y(0, uj) + w{0, n^, oo))] + P^G(ft,w, y{0, n^, uj)). 

{-A - z{9rUj))P''yiO, Mc(r, v, uj),uj) + -^P"y(0, u,{r, v, uj),uj) 

ar 

= {-A - z{9rUj))P''y{0, u,{r, v, uj),uj) + P"y„(0, u,{r, v, uj),uj) 

[AUc + z{6rUj)Uc + P^GiOrUJ, 1/(0, Oj) + w(0, Uc, Oj))] 

+ P"i/r(0,nc(r, 
= M(P"y) - P"y,(0, Uc{s, v, u),uj) [P^G(a;, y{0, m„ u)) 

- P''G{9rdrOO, 1/(0, Mc, + w(0, Mc, t^))] + P^'GiOrUJ, 1/(0, lie, w))- 

From the above calculations 



7^^(1/(0, i;, w) + ii;(0, l^ w), 0) - 1/(0, t;, w) 



.0 




POO 




/ exp 

/ — oo 


-Ar + 1 z{9iu) dl 

•J r 


Q" dr + / exp 
Jo 


-Ar+ j z{9iu)dl 



where and is 



Q'{y,w,u) = P'y^{0,Uc{r,v,u),uj)[P''G{9rUJ,y{0,Uc,u)) 

- P^GiOrU, 1/(0, n„ u) + if;(0, u^, cu))] - M{P'y) 

+ [P'GiOrUj, 1/(0, n„ a;) + ii;(0, ix„ a;)) - P^G(ft,a;, i/(0, n,, a;))], 

Q"(i/, 1/;, oo) = M{P-y) - P"i/,(0, ii,(r, i;, a;), [P^G'(e,a;, i/(0, iz„ a;)) 

P^GiOrto, 1/(0, He, u;) + ii'(0, Ue, 0;))] - [P"G(^r-c^, 1/(0, i/e, 0;) 
+ il;(0, Ue, w)) - P"G(ft,a;, i/(0, n,, w))]. 

Since gv{0, u) = and (/ is G^ in a neighborhood of origin in iP^, for a small 
random variable e{u), 

\gi,{v,uj)\ < e{u) 
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if small enough. From M{y){0,Uc,uj) = Ci{uj)\uc\'^ for a positive ran- 
dom variable Ci{uj) and all u^. G H'^, we obtain \Q^\ < C2{uj)\uc\'^ and 
\Q^\ < C3{u)\uc\'^ for some random variables 6*2(0;), Cs{uj),. Since Uc is 
the solution of equation (35), from the Gronwall's inequality we conclude 
\uc\ < 6*4(0;) exp(— cr) |w| for a positive constant c and random variable C*4(cj). 
Since the random variables 

exp( / z{6iuj) dl) dr < +00, / exp( / z{6icu) dl) dr < +00, 

-00 J r Jo J r 

SO |5w;(0,w,a;)| = | J"(w;(0, ^^ w) + ?/(0, w), 0) - ?/(0, a;)| < 6*5(0;) lu^l'? < 
K{uj)\v\^ by selecting suitable random variables 6*5(0;), A'(o;). □ 

7 Examples show our results 

Example 7.1 shows that we can use computer algebra [29] to directly compute 
stochastic center manifolds in a wide class of spdes. Example 7.2 establishes 
that a semilinear damped wave equation has a stochastic center manifold. 
Lastly, Example 7.3 explores where the stochastic center manifold 

itself is infinite dimensional. 

7.1 Stochastic reaction-diffusion modelling 

Consider the stochastic reaction-diffusion equation 

Ut = Uxx + u — au^ + auoW with u(0, t) = u{7r, t) = , t > . (36) 

We set the space H := L^(0, vr), the operator Au = d'^u/dx'^ + u with domain 
D{A) := H'^{0, vr) fl -f^o(0, tt). Then the spectrum of the operator A is 

a{A) = {0, -3, -8, -15, . . .} = {1 - : n = 1, 2, . . .}. 

The corresponding complete and orthogonal eigenfunctions are sin nx , n = 
1,2,.... The operator A satisfies the exponential trichotomy (Definition 1): 
the center subspace = span{sinx}, the stable subspace H'^ = span{sinnx, 
n = 2,3, . . .}, and the unstable subspace if" = 0. 

Recall that the random variable z{6tuj) is the stationary solution of the 
Ornstein-Uhlenbeck SDE (6). In order to transform the random parabolic 
equation (36), from the coordinate transform (7) we let the transform u = 
u* exp{az{dtuj)). For simplicity omit the * from here on, the transformed 
random parabolic equation is 

Ut = Uxx + u + az{9tuj)u — aexp[2az{9tuj)]u'^. (37) 
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Select the nonlinearity parameter a small enough such the conditions of 
Theorem 5 is satisfied. Then Theorem 5 assures us that there exists a slow 
(center) manifold 

{oo 
s sinx + c„(s, u) sinnx 
n=2 

where Cn{s, co) = 0(^s^) as s — > for to E . The random dynamical system 
on the slow manifold M'^[u) is governed by the following equation. 

ds 3 

— sinx = saz{9tu) sinx — -exp[2az{9tCo)]as^ sinx + 0(^s^). (38) 

Equation (38) has a common factor of sin x in each term, it follows that 
ds 3 

— = saz{etu) - -exp[2az{9tu)]as^ + 0{s^). 

We now calculate the slow manifold h'^{s sinx, u). For u E fl, Cn{s, u) is odd 
functions of s since —u is a solution of equation (37). The coefficient Cn{s, u) 
is zero for odd term, so 

h'^{s sinx, u) = 03(3, co) sin3x + 0[s^) . 

From C3{s,uj) = we obtain 

aexTp[2az(9tUj)]{s sinx + C3 sin 3x)'^ = as^ exp[2az(9tuj)] sin'^ x + 0[s^) . 
So 

3 

P''G{9tUJ, s sinx + h") = -as^ exp[2az{9tuj)]smx + 0{s^) , 

P'G{9tU, s sin X + h"") = -^as^ exp[2az{9tu)] sin 3x + 0{s^). 

From the formulation of approximation of center manifold, 

M{h'){v,uj) = h'^{v,uj) + hl{v,uj)[P^Av + z{9tuj)v + P''G{9tio,v + h')] - Ah^ 

- z{9tuj)h^ - P'G{9tu;, v + h"") - P''G{9tu, v + h"), 

there exists an stochastic approximation of center manifold of equation (37). 

Now we directly use computer algebra to compute the stochastic center 
manifold [29]. Define the Ornstein-Uhlenbeck processes (f)k{t) = exp[— (/c^ — 
l)t] ^iy(t), then 

u = s sinx + —as'^ sin 3x — — as^(703(t) sin3x + 0(3^,0^, cx^) (39) 
32 16 

ds ■ 3 

such that — = saW(t) - -as^ + 0(s^,a^,a^). 

dt 4 V ' ' / 
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7.2 Semilinear damped wave dynamics 



Wang and Duan [36] discussed a semilinear wave equation with a two dimen- 
sional center-unstable manifold. We establish that the semilinear damped 
wave equation has a one dimensional center manifold. 
Consider the spde 

utt + ut = ^u^^ + u + f{u) + uoW, xG(0,27r), t>0, (40) 

with u{0,t) = u{2TT,t) = 0. W{t) is the standard M-valued Wiener process 
on a probability space {Q,J^,F). The nonlinearity term / is smooth and 
/(O) = 0. 

Let the space H = Hq{0,2it) x L^(0, 27r). Rewrite the spde (40) as the 
following first order stochastic evolution equations in H: 



du = V dt , 

dv = [^u^^ + u - V + f{u)]dt + uo dW(t). 



(41) 
(42) 



(43) 



First we prove that the system (41)-(42) generates a continuous random 
dynamical system in H. Let "^lit) = u{t), '^2it) = v(t) — u{t)z{9tuj), where 
z{u) is the stationary solution of (6) with /i = 1. Then equations (41) and 
(42) are transformed to the following random evolution equations 

rf^^ = + ^iz{0tuj)] dt , 

d^2= [^^ + ^i-^2 + f{'^i)]dt+[z{etUj)^,-; 

(44) 

Let ^(t,w) = (^i(t, w), ^2(t, w)) e H , then Lp{t,'^{0,uj),uj) = ^(t,w) de- 
fines a continuous random dynamical system in H. Notice that the stochastic 
system (41)-(42) is conjugated to the random system (43)-(44) by the home- 
omorphism 

T{u, {u,v)) = {u,v + uz{u)), {u^v) E H 

with inverse 

T~'^{ijj, (u, v)) = {u,v — uz(uj)), (u, v) E H . 

From Lemma 3, (p(t,u, {uo,vo)) = T{9tU,ip{t,u,T~^{u, {uo,vo)))) is a ran- 
dom dynamical system generated by (41)-(42). Define 



A = 



0, 



4 5x2 



id ■ 
-id 



+ id, - 

ziOtoo), 

-z^etco), ~z{etu) 
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where id is the identity operator on the Hilbert space L^(0,27r). Then the 
random system (43)-(44) can be written as 

— = + Z{9tuj)^ + F{^). (45) 

(Jjlj 

The operator A has the eigenvalues 

and corresponding eigenvectors are = {1, 5^) sinkx . The operator A 
has one zero eigenvalue 5^ = , one positive eigenvalue 5]^ = | and the 
others are all complex numbers with negative real part. By using a new 
inner product [3(3], we define an equivalent norm on H with the eigenvectors 
{(1, 5^) sinkx} orthogonal in this norm. 

A calculation yields that, in terms of the new norm, the semigroup S{t) 
generated by A satisfies the exponent trichotomy with a = ^ + e,'j = ^ + e 
/3 = I — e for parameter e small enough, and dim = dim = 1 . 

Let \1/ = rqf + sq2 +y , y G P**if . For application the theory of stochastic 
center manifolds, projecting the random system (45) to the center subspace 
span{g^} by the projection P'^, 

sqt = Z{etuj)sq+ + P'F^OtUJ, rqf + sqt + y). (46) 

By a cut-off function. Lip F is sufficient small. From Theorem 5 we obtain a 
center manifold Miu) which we denote as the graph of a random Lipschitz 
map h^{sq2,uj). 

The asymptotic behavior of the random system (45) is governed by the 
one dimensional equation 

sg+ = Z{etuj)sq+ + P^F^OtU, sq+ + h%sq+, oo)). (47) 

From Theorem 10, we know that the long time asymptotic behavior of the 
random system (45) is the same as the equation (47). 

7.3 An infinite dimensional center manifold 

There are very few rigorous results involving infinite dimensional center man- 
ifolds except the existence theorem of Gallay [20] and the cylindrical analysis 
of Mielke [23]. However, applied mathematicians and engineers implicitly, in 
effect, assume an infinite dimensional center manifold whenever they invoke 



29 



the method of multiple scales for slowly varying solutions [?, ?, e.g.]. As a 
basic example let's consider the coupled spdes 

dzi • dv d'^v 

= au - uv + auoW , ^ = -v + — ^ + - 2}Ca(u^v) + avoW , (48) 
ot ot ox^ 

with Neumann boundary conditions = at 2; = 0, vr , and the operator 
/Ca = (1 + 2a — dzz)~^ ■ For example, for the specific case a = 0, 



JCqV = / Kq{z, ()v{C) d( where -f^o 



cos(7r— (^) cos z 

sin-TT 
cos(7r— 2:) cos f 



, z < C. 

Sm TT ' ^ 



sm-TT 



z>C, 



for which esc vr < Kq < cot tt and Kq d( = 1 . 

Let the transform u = u* exp((T2;(6'ta;)), where the random variable z{6tU}) 
is the stationary solution of the Ornstein-Uhlenbeck SDE (49). 

dz + zdt = dW. (49) 

Transform the SPDEs (48) to the following system of spdes with random 
coefficients: 

Qu 

— = au + az(6tui)u — exp(a z(6tui))uv , 
ot 

dv d'^v 

— = -w + — + az{6tuj)v + ex.Y){az{0tUj))u'^ - 2 ex.Y){2az{0tUj))lCa{u'^v) , 

(50) 

If a = 0, from Theorem 5, there exists an infinite dimension stochastic 
slow manifold. Assume the stochastic slow manifold is f = h'^{u,uj), then 
from Theorem 12, the exact slow manifold satisfies 

h^u, u) + hl^iu, u) [az{9tuj)u — exp{az{u))uh'^{u, u)] 

+ (h'^(u, co) — ——'^(u, Lo)) — az{6tuj)h^{u, uj) — exp(az(9tco))u^ 

OX'' 

+ 2exp{2az{etLj))]Ca{u^h%u,uj)) = 0, 

for cr small enough. From the construction of operator JCq, the above system 
of SPDEs (50) has a stochastic slow manifold of 

V = h{u,u) = exp{a z{9tUj))ICou'^ . 

By the equation (28), the evolution on the center manifold is governed by 
the SPDE 

du 

— = az{9tuj)u — exp{2az{9tUj))u]Cou'^ . 
ot 
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